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ABSTRACT 

Future large arrays of telescopes, used as intensity interferometers, can be used to 
image the surfaces of stars with unprecedented angular resolution. Fast-rotating, 
hot stars are particularly attractive targets for intensity interferometry since shorter 
(blue) wavelength observations do not pose additional challenges. Starting from real¬ 
istic surface brightness simulations of fast-rotating stars, we discuss the capabilities 
of future intensity interferometers for imaging effects such as gravity darkening and 
rotational deformation. We find that two-telescope intensity correlation data allow 
reasonably good imaging of these phenomena, but can be improved with additional 
higher order (e.g. three-telescope) correlation data, which contain some Fourier phase 
information. 


Key words: instrumentation: interferometers-techniquesdhgh angular resolution- 
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1 INTRODUCTION 

Over half a century after t he success of the Narr abri Stel¬ 
lar Intensity Interferometer (iHanburv Brownl[l974T ). a future 
facility: the Cherenkov Telescope Array (CTA), has been 
proposed to be used in part as an Intensity Interferometer 
( Le Bohec fe Holderll2006l : iDravins et al.ll2013l ; INunez et al.l 
l2012al l. When used as an optical Intensity Interferometer, 
CTA will provide thousands of simultaneous baselines as 
large as a kilometer, corresponding to sub-milliarcsecond 
angular resolution at visible wavelengths, and thus enable 
to resolve the surfaces of many nearby stars. Since intensity 
interferometry is insensitive to atmospheric turbulence and 
optical imperfections, observing at short blue wavelengths, 
where Cherenkov telescopes are particularly sensitive, does 
not impose additional challenges. Also, Intensity Interferom¬ 
etry has been shown to be particularly sensitive to hot stars, 
which emit much of their radiation in the blue. Such observa¬ 
tions would complement amplitude interferometry observa¬ 
tions, which are typically done at longer wavelengths. Some 
numerical simulation studies of future intensit y interferom¬ 
eters s how that stellar surfaces can be imaged (INunez et al.l 
2012a||b|). and recent laboratory e xperiments performed by 
Dravins, Lagadec fe Nufie5 (!2015l l have resulted in the first 
optical images from intensity interferometry data. 

In this simulation paper we investigate the capabilities 
of CTA used for Stellar Intensity Interferometry (SII) obser¬ 


vations of hot, fast-rotating stars (spectral types O, B, and 
A) at blue wavelengths. These stars are those presenting the 
highest rotation rates in the H-R diagram for non-degenerate 
objects. The most apparent effects of fast-rotation on the 
stellar surface are (1) the deformation of the outer layers of 
the star (rotation flattening), and (2) the latitudinal depen¬ 
dent emitted flux (gravity darkening), both consequences 
of the balance between significant centrifugal and gravita¬ 
tional forces within the star. Depending on the rotation rate, 
these two effects can lead to significant modifications in the 
specific intensity distribution of the star as seen by a dis¬ 
tant observer. SII will allow performing morphological stud¬ 
ies studies, e.g. of rotational flattening, as well as mapping 
the bri ghtness (temperat ure) distribution on the stellar sur¬ 
face (INunez et al.l l2012bh . The effect of stellar rotation on 
the surface brightness distribution affects our understand¬ 
ing of radiation transport, which may depart from simple 
diffusion, and affect the way elements can be re-distributed 
within the star. Therefore, imaging at high angular reso¬ 
lution may in turn furt her our understan ding of the time 
evolution of these stars dDelaa et al.ll20ljh . 

The outline of this paper is the following: In Section 
[2] we introduce the model used to simulate the surface ra¬ 
diance distribution of fast rotating stars. In Section [3] we 
discuss two- and three-point intensity correlations as well 
as their Signal-to-Noise-Ratio (SNR). In Section [4] we pro¬ 
vide a detailed description of the data simulation. Simulated 
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and laboratory image reconstructions have so far relied ex¬ 
clusively on two-point intensity correlation data, and one of 
the main goals here is to investigate the benefit of addition¬ 
ally using three-point correlations. In Section [o] we describe 
our data analysis approach, focusing on model-independent 
image reconstruction from the simulated data, and present 
our results using two- and three-point intensity correlation 
(simulated) data. We conclude in Section [ 6 ] 


2 FAST-ROTATING STARS 


The study of fast-rotating stars has largely benefited from 
the high angular resoluti on capabilities of modern ampl i¬ 
tude interferometry (e.g. Domician o de Souza et al. 1 12014 
iMonnier et alJ 120121 : Ivan Bell? 20121 '). These works showed 
that the photospheric structure of rapidly rotating, non¬ 
degenerate, single stars of intermediate to high masses are 
well described by a model based on the following assump¬ 
tions: solid rotation (angular velocity fl), mass M concen¬ 
trated in the center of the star, and gravity darkening. From 
these assumptions, the stellar surface of the star is deformed 
(flattened) by rotation and follows the Roche equipotential 
(gravitational plus centrifugal), 


* (*) = - 


GM fl 2 R 2 (9) sin 2 6 GM 
R(0) 2 ~ 


2 


(i) 


where 9 is the co-latitude, f? eq and u eq (= f 2 f? eq ) are the 
equatorial radius and rotation velocity. Solving the Roche 
equipotential provides the co-latitude-dependent stellar ra¬ 
dius R(9). The local effective gravity is obtained directly 
from the gradient of Roche equipotential , i.e., 


ffe ff(0) = |-Vtt(0)| • ( 2 ) 

Gravity darkening is included by relating g e s(6) to the 
local effective temperature T e ff($) in a generalized f orm o f 
the gravity-darkening law proposed by Ivon Zeipell dl924f) . 
i.e., 

T eff (0) = Kg p eS (9) , (3) 


where A is a proportionality constant, and /3 is the 
gravity-darkening coefficient (/? = 0.25 in the origi¬ 

nal von Zeipel gravity darkening). Deviations from j3 = 
0.25 are interesting since they are indic ative of con vec- 
tive atmospheric physic s (e.g. [5 = 0.08 (jLuc 5I1967D, or 
something in between (IZhao et al.| |20 09l; IChe et al.l l 201 ll ; 


lEspinosa Lara fe Rieutordf 2011 : Clarelll2015l ). This model. 

hereafter called RVZ (for Roche-von Zeipel) model is used in 
the next sections to simulate SII observations of fast-rotating 
stars. 


3 INTENSITY INTERFEROMETRY 

In a Stellar Intensity Interferometer, the star-light received 
at each aperture is converted into a time-dependent intensity 
signal with fast electronics. Time averaged correlations are 
found between telescope pairs (and possibly triplets, etc.) 
with analog correlators or off-line digital signal processing. 
Since the tolerance on optical imperfections/delays is set by 
the electronic time resolution, Intensity Interferometry can 


be insensitive to optical delays of the order of a meter when 
using detectors with a response time of a few nano-seconds. 

These intensity correlations in turn provide information 
about the degree of coherence of light, i.e. the interferometric 
visibility. The time-averaged two-point intensity correlations 
( IiJj) are related to the interferometric (complex) visibility 
7 ij as 


— 1 + \lij\ 


( 4 ) 


(SIA 
(Ii > (Ij) 

Since the phase of the complex visibility is not mea¬ 
sured with two-point intensity correlations, image recon¬ 
struction from two-point correlations is non-trivial and has 
benefited from the use of several p hase retrieval methods and 
image reconstruction techniques (|Hohnes_fc Belenikiil 120041 : 
iNufiez et al.ll20f~2al : IStrekalov, Kulikov fe YiJ 2014 1. 

However, it can be shown that higher order correla¬ 
tions can retain some Fourier phase information. In the case 
of three-point-correlations, it can be shown that (see e.g. 
dMalvimat, Wucknitz fe Sahall2014l l 

/r \ //wl v = l + l7ij| 2 + l7j'fc| 2 + l7fc*| 2 +2Re[7ij7jfc7fci] .(5) 
VR Vi) Vk) 

The phase of the triple product in last term of this 
equation is known as the closure phase, i.e. 


(j> c = Arg [jij'TjkJki] ■ ( 6 ) 

The closure phase (rf> c ) has the property of being insen¬ 
sitive to the atmospheric phase since baselines ij, jk and ki 
form a closed loop. Strictly speaking, three-point intensity 
correlations only allow measuring the cosine of the closure 
phase. In Section [5TTI we discuss strategies for reconstruct¬ 
ing images using two- and three-point intensity correlation 
measurements. While some simulated studies on image re¬ 
construction have been performed using two-point intensity 
correlations, closure-phase intensity-correlation imaging is 
still in embryonic form. 


3.1 Signal-to-Noise 


The signal-to-noise (SNR( 2 )) for a two-point intensity cor¬ 
relation measurement depends on the degree of coherence 7 , 
the area A of each of the light receivers, the spectral density 
n (number of photons per unit area per unit time, per fre¬ 
quency), the quantum efficiency a, the electronic bandwidth 
A/, and the observat ion time t. The explicit SNR expression 
can be shown to be dHanburv Brownlll974T ) 

SNR( 2 ) = n(\, T,m v )Aa | 7| 2 \JAft/2. (7) 


Note that SNR ( 2 ) is independent of the spectral band¬ 
width, imposing no additional difficulty of observing at in¬ 
dividual spectral lines. A more deta iled discussion of th e 
sensitivity of CTA for SII is given bv I Nunez et al.l (|2012al l. 
and according to their simulations, tens of hours of integra¬ 
tion time are required for observing bright ( m v ~ 6 ) stars 
at blue wavelengths. 

A study of the SNR( 3 ) fo r three-point correlations has 
been recently performed by iMalvimat. Wucknitz fe Sahal 
(l2014ll and further investigated by Wentz fe Sahal ( 2015lh 
who find 


SNR ( 3 ) = | 7 ij 7 jfe 7 fci|[n(A, T, m„)A] 3 / 2 A/ v / t/Au, 


(8) 
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where 7 jj is the degree of coherence for the baseline i — j, n 
is the spectral density, A is the individual telescope area, A/ 
is the electronic bandwidth, t is the observation time, and 
Ai/ is the optical bandwidth cAA/A 2 . Equations [7] and [ 8 ] will 
be used for adding noise to the simulated data presented in 
Section 0] 

If we assume individual CTA telescope a reas of ~ 43 m 2 
(rough ly ~ 3.7 m dishes consistent with iBernlohr et al.l 
(I2013fl 'l. an electronic bandwidth of ~ 2 GHz, a quantum ef¬ 
ficiency of ~ 60%, A = 545 nm, and an optical bandwidth of 
~ 10 11 s ' 1 for 50hrs of observing time on a 3 rd magnitude 
star (corresponding to n ~ 10 ~ 5 m~ 2 s ~ 1 Hz ^ 1 , depending 
on stellar temperature), we obtain SNR( 2 ) ~ 2000 and a 
SNR ( 3 ) ~ 3. Both SNR( 2 ) and SNR ( 3 ) can benefit from 
using multiple spectral channels since the SNR improves as 
the square root of the number of spectral channels. In Sec. 
I5.2l we estimate the number of spectral channnels needed to 
benefit from the additional use of three-point correlations. 


4 DATA SIMULATION 

4.1 Simulation of the Pristine Images 

The simulation process starts by generating an angular sur¬ 
face brightness distribution (or intensity maps) in a certain 
wavelength range (see Figure [T]). 

These images are generated with the numerical model 
CHARRON (Code for High Angular Resolution of Rotating 
Objects in Nature) that we shor tly describe below. A more 
detail e d des cription is given by iDomiciano de Souza et alJ 
li 20021 . 2012). The adopted physical model for the fast- 
rotating star is the RVZ model summarized in Sect. [2] 

CHARRON is an IDL-based numerical model, oriented 
to deal with spectro-interferometric data, which means that, 
from the input physical parameters, it computes monochro¬ 
matic intensity maps of the stellar photosphere. From 
the Fourier transform (FT) of these maps one can ob¬ 
tain spectro-inferferometric observables (e.g. visibility am¬ 
plitudes, differential phases, closure phases). Numerically, 
the stellar photosphere (Roche potential surface) is divided 
in a predefined grid with nearly identical surface area ele¬ 
ments (typically ~ 50 000). From the RVZ physical equa¬ 
tions it is possible to assign different physical quantities 
to each surface grid element j-. radius Rj, rotation velocity 
Vj (= QRj), effective gravity g e s,j, and effective temperature 
T e s,j ■ 

A local specific intensity from a plane-parallel atmo¬ 
sphere Ij(= Ij(g e fc : j, T e ff y j, A(u P roj,j), Uj) is then associated 
to each surface element. The wavelength A(uprojj), corre¬ 
sponding to the local specific intensity, is Doppler-shifted to 
the local rotation velocity projected onto the observer’s di¬ 
rection Uproj ,j ■ The parameter gj is the cosine between the 
normal to the surface grid element and the line-of-sight (limb 
darkening is thus automatically included in the model). Both 
Vprojj and gj depend on the chosen inclination i of the ro¬ 
tation axis. 

The local Ij can be given by an analytic equation (in¬ 
cluding for example continuum flux, line profiles, limb dark¬ 
ening) or as an output of stellar atmosphere models calcu¬ 
lated with radiative transfer codes. The results presented 
in this work are mainly based on Ij interpolated on a pre¬ 
calculated grid of specific intensities. These pre-calculated 


Table 1 . Stellar parameters adopted to create the model images 
shown in Fig. [I] The adopted parameter values roughly corre¬ 
spond to a fast-rotating B type star. 


Model parameter 

Value 

Stellar mass M 

6 Mq 

Inclination angle of rotation axis i 

0°,45°,60°,90° 

Equatorial radius R eq 

11 Rq 

Polar radius R p 

7.83 Rq 

Equatorial rotation velocity u eq 

290 kms -1 

Gravity-darkening coefficient /3 

0.20 

Surface averaged mean temperature T e ff 

15 000 K 

Polar eff. temperature T p 

17 819 K 

Equatorial eff. temperature T eq 

11180 K 



1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 

Surface Brightness in Arbitrary Units 


Figure 1. Simulated Pristine images at 545 nm. Each image cor¬ 
responds to a different inclination angle of the stellar rotation 
axis, where the angle is taken with respect to the line of sight. 
These models of fast rotators were calculated with the CHAR¬ 
RON code summarized in Sect. [4] 


specific intensities were obtained from th e spectral synthe¬ 
sis code SYNSPEC dHubenv fc Lani20ljh and the ATLAS9 
stellar atmosphere models ( Kuruc j 1979ft . The total inten¬ 
sity map of the star I\ is obtained from the juxtaposition 
of each Ij for the apparent surface grid elements. Examples 
of these intensity maps (our pristine images) are given in 
Figure [l] 


4.2 Simulation of Intensity Interferometry data 

In order to simulate SII data, we follow the procedure de¬ 
scribed in detail by I Nunez et al.l (1201281 1. First we set the 
spatial scale of the image: we concentrate on the case of a 
bright star of m v = 3, which can have an angular diame¬ 
ter of ~ 0.5 mas as long as the temperature of the star is 
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Figure 2. Telescope array configuration ( jBernlohr et alJl2013h 
used for the simulations. 
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more thaifi ~ 10, 000 K. Larger angular diameters than our 
choice of 0.5 mas are acceptable as long as they are not fully 
resolved by baselines smaller than ~ 100 m. If needed, the 
pupil of the largest telescopes (~ 24 m, presumably located 
towards the center of the array) could be segmented into sub¬ 
pupils to allow observations with angular diameters much 
larger than 0.5 mas. Then we take the Fast Fourier Trans¬ 
form (FFT) of the pristine image, and sample the data in 
the ( u , v ) plane as determined by the telescope array con¬ 
figuration and the o bserving wavelength . A possible CTA 
array configuration (jBernl ohr ct. al .1120131 ) is shown in Fig¬ 
ure [H consisting of 74 telescopes of 4 m and 7 m telescopes, 
spread out over 2 kilometers. Such an array provides 2701 
baselines pairs, ranging from 14 m to 2 km, and 70300 tele¬ 
scope triplets. Two-point and three-point correlation mea¬ 
surements are simulated by sampling an FFT over the base¬ 
lines (pairs) and triplets provided by the array, using equa¬ 
tions [4] and [5] Gaussian noise is introduced according to 
equations [3 and [8] 

We simulate intensity correlation data at 545 nm for 
each of the pristine images shown in Figurc[T] For simplicity 
we assume that the star is observed at zenith for a total of 
50 hrs, and that baselines do not drift due to earth’s rotation, 
achieved by performing short nightly observations of, say, a 
few minutes. In Fig. [3] we show projections of the two-point 
intensity correlation data corresponding to m v = 3, using 
a single spectral channel. The data simulation parameters 
are summarized in Table [2] To simulate data with multi¬ 
ple spectral channels close to each other, we simply reduce 


1 To estimate the temperature corresponding to a particular 
m v and angular diameter, we use the fact that m v — mo = 
—2.5 log{F/F q), where Fq is the flux corresponding to a star of 
mj and F is the flux of our simulated star. Last, we recall that 
the flux is proportional to 9 2 T 4 , where 9 is the angular diameter 
and T is the effective temperature. 


Table 2. Data simulation parameters. 


Data Simulation Parameter 

Value 

Equatorial Angular Diameter 

0.5 mas 

Apparent Visual Magnitude 

3 

Observation Time 

50 hrs 

Wavelength 

545 nm 

Optical Bandwidth per spectral channel 

10 11 s” 1 


the error bars by the square root of the number of spectral 
channels. 


5 DATA ANALYSIS AND INTERPRETATION 

As in amplitude interferometry, the two-point intensity cor¬ 
relations provide physical information on fast-rotators, such 
as equatorial angular diameter and rotation flattening. As 
shown in Figure [3] the high spatial frequencies attained by 
SII allow to probe the second and third visibility lobes, 
which are more sensitive to limb and gravity darkening than 
the first lobe alone. The three-point correlation data can be 
sensitive to asymmetries in the stellar brightness distribu¬ 
tion, mainly due to gravity darkening in the case of fast 
rotating stars. 

A commonly adopted strategy to extract the physical 
information from the data is to apply a model fitting pro¬ 
cedure that allows to estimate some model parameters and 
their uncertainties. For modern amplitude interferometers 
the uncertainties on the model parameters range typically 
of t he order of ~ 1% — 10%, depending on the paramete r 
(e.g. iDomiciano de Souza et all 120141 : iMonnier et al.ll2012h . 
Since the (u, v) coverage of CTA can be so dense, we rather 
adopt model-independent image reconstruction strategies 
described in Sec. EQ An important point we will investi¬ 
gate is the imaging benefit of using three-point correlations 
in addition to two-point correlation data. 


5.1 Image reconstruction strategy 

We test model-independent image reconstruction tools to 
analyze the data displayed in Fi gure [3] We first follo w the 
strategy described in detail by (I Nunez et al.l l2012bl f : The 
first step consist in estimating the phase from Fourier magni¬ 
tude data alone, using the “Cauchy-Riemann” phase recon- 
struction method, based o n the theory of analy tic functions 
dHolmes fe Belen kiil 120041 : INunez et al.l l2012al l. The main 
idea behind the Cauchy-Riemann method is that the Fourier 
Transform of an object of finite extent is an analytic func¬ 
tion. Therefore, one can use the Cauchy-Riemann equations 
(in polar coordinates) to relate the Fourier magnitude to 
the phase. Since the Cauchy-Riemann algorithm relies on 
the calculation of derivativ es, an analyt i c funct ion is fitted 
to the data as described in INunez et al.l d2012a| j. 

The first estimate can be further improved by using 
two different algo rithms: the “Gerchberg-Saxton” algorithm 
dGerchbereHl972T ). which consists in going “back and forth” 
between image and Fourier space, each time imposing gen¬ 
eral restrictions in each domain. In our case, the “Gerchberg- 
Saxton” restrictions in image space consist in setting setting 
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i=0 


i=45 



Figure 3. Two-point intensity correlation data at 545 nm for m v = 3 and 50hrs of observation time for the different inclination angles of 
the fast-rotating star. The data labeled as 0 deg corresponds to baselines parallel to the x-axis in Fig. [2] the other data points correspond 
to baselines at 45° and 90° with respect to the x-axis in Fig. [7] The data are binned every 10 m for visual clarity. 


pixels outside a circular region to zero, where the size of the 
circular region is determined by the smallest closest zero 
of the Fourier Magnitude. The restriction in Fourier space 
is that the Modulus of the Fourier transform of the recon¬ 
structed image should be replaced by the analytic function 
fit. We note that a poor fit can already introduce artifacts 
in the reconstructed image. 

The reconstructed image can be further improved with 
yet another algorithm: the “M ulti-aperture Im age Recon¬ 
struction Algorithm” (MiRA) dThiibautl 120091 ') : an opti¬ 
mization algorithm which maximizes agreement between the 
Fourier transform of the reconstructed image and the data. 
The MiRA is particularly good for removing artifacts that 
may have arisen from the initial fit of the data. The reason 
for using a seemingly complex chain of algorithms, and not 
a single algorithm, has been justi f ied em pirically via simu¬ 
lations performed bv lNunez et al.l (l2012bi l. and the different 
algorithms seem to complement each-other. In the context 
of this paper, the MiRA is also particularly useful for us¬ 
ing the three-point correlation (closure-phase) data, and to 
test if such higher-order correlation measurements are in¬ 
deed useful for image reconstruction. 

Recall from eq. [5] that a three-point intensity correla¬ 


tion allows measuring the cosine of the triple product of 
the degrees of coherence of in a closed loop of baselines. 
The phase of the triple product (eq. [6J is the closure phase 
known in amplitude interferometry, but only the cosine of 
the closure phase is accessible with intensity interferometry 
measurements. This results in having a sign (and modulo 
27 t) ambiguity in the closure phase measurement. Provided 
sufficient SNR, the sign ambiguity should be removable by 
performing observations at two (or more) close wavelengths 
separated by AA, and estimating the change in the three- 
point correlation for a given AA. Measuring three point cor¬ 
relations at two close wavelengths yields two closure phase 
equations of the form of Equation [S] If the object’s surface 
brightness distribution does not differ largely between the 
two different wavelengths, the two closure phase equations 
should mainly differ by an additive term proportional to the 
ratio of the two wavelengths (in a first order approximation). 
This should allow removing the sign ambiguity by choosing 
the sign of the closure phase from one of the two closure 
phase equations, and then checking if it is consistent with 
the other equation. 

The problem of removing the sign ambiguity in the clo¬ 
sure phase and how it affects the SNR deserves further in- 
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Table 3. Mean Square Error quantifying the difference between 
the pristine images (left column) and the reconstructed images 
(top row) without closure phase data , using a single spectral chan¬ 
nel. The images used for computing the Mean Square Error are 
shown in Figures [l] and [4] The minimum MSE of each column is 
shown in bold-face. Note that the reconstruction of the case of 
i = 45 can be confused with a pristine image with i = 0. 


Prist\Rec 

i=0 

i=45 

i=60 

i=90 

i=0 

1.37 X 10~ 3 

1.24 X icr 3 

5.50 x 10~ 3 

5.39 x 10“ 3 

i=45 

6.44 x icr 3 

1.43 X 10~ 3 

1.43 x 10~ 3 

2.00 x 10 -3 

i=60 

1.24 X icr 2 

5.25 X 10 -3 

8.43 X 10~ 4 

2.40 x 10 -3 

i=90 

9.92 x icr 3 

3.36 X 10“ 3 

1.92 x 10~ 3 

9.20 X 10 -4 


vestigation and will be treated in a future paper. In the 
simulations presented here, we assume that the sign ambi¬ 
guity has already been removed with the approach described 
above, and we have not included this phase degeneracy re¬ 
moval in our simulation^. The main interest is to see how 
three-point correlations with (much) lower SNR than two- 
point correlations can improve image reconstructions. 


5.2 Image reconstruction results 

In this section we compare image reconstructions using two- 
point correlation data with reconstructions using additional 
three-point correlation data. We first investigate image re¬ 
construction quality when using only two-point correlation 
data for a single spectral channel, and find reasonably good 
image reconstructions as shown in Figure [4] which may al¬ 
ready allow to infer the rotational flattening and the inclina¬ 
tion angle in some cases. To quantify image reconstruction 
quality, we use the Mean Square Error (MSE) between the 
pristine image and the reconstructed image. In order to be 
more certain that we can correctly determine the stellar pa¬ 
rameters from an image reconstruction, we compare each 
reconstruction in Figure [4] with all the different pristine im¬ 
ages in FigurcQ] and provide the values of the MSE between 
each reconstructed images and pristine images in Table [3] 
For each reconstruction, the lowest MSE should be found 
for its corresponding pristine image. As can be seen in the 
Table, the reconstruction for i = 45° is more similar to the 
pristine image with i = 0°, i.e. we have reconstructed at 
least one image for which we cannot determine the correct 
inclination angle and gravity darkening. 

Next we investigate whether or not three-point correla¬ 
tion data (in addition to two-point correlations) over a single 
spectral channel improve image reconstructions. We showed 
in Sec. o that SNR( 3 -) ~ 3 for a single spectral channel, 
so in principle, we should be able to see improvements when 
using three-point correlation data. We also expected image 
reconstructions to improve due to the much larger number 
of telescope triplets than pairs, but we do not find a sig¬ 
nificant improvement using the analysis method presented 
in this paper, at least for the case of this object observed 
with this array configuration. Clearly, there is still room for 
improvement on the data analysis method. 


2 All the codes used for our simulations are available upon re¬ 
quest. 
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Surface Brightness in Arbitrary Units 

Figure 4. Image reconstructions without three-point correlation 
information. 
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Figure 5. Image reconstructions using three-point correlation 
information. 


MSE for the case of i=45° 
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SNR( 3 ) Enhancement Factor 

Figure 6. Mean Square Error for the image reconstruction cor¬ 
responding to the case i = 45° with its corresponding pristine 
image as a function of the enhancement factor for SNR^y The 
horizontal line is the lowest MSE found in the second column of 
Table [3] A reliable improvement in the image reconstruction can 
be seen for an enhancement factor approximately greater than 
- 30. 
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Table 4. Mean Square Error between the pristine images (left col¬ 
umn) and the reconstructed images (top row) with closure phase 
data, using 1000 spectral channels. The images used for comput¬ 
ing the Mean Square Error are shown in Figures [l] and [5] The 
minimum MSE of each column is shown in bold-face. Note that 
it is unlikely to confuse a reconstruction with the incorrect pris¬ 
tine image. 


Prist \Rec 

i=0 

i=45 

i=60 

i=90 

i=0 

1.17 x 10~ 3 

1.44 X 10“ 3 

4.62 x 10“ 3 

5.67 x 10" 3 

i=45 

2.32 X 10 -3 

1.01 X 10~ 3 

4.47 x 10“ 4 

2.32 x 10 -3 

i=60 

1.15 X 10 -2 

4.27 X 10 -3 

7.16 X 10~ 4 

2.97 x 10 -3 

i=90 

9.24 x 10 -3 

3.16 x 10 -3 

1.60 x 10 -3 

1.02 X 10“ 3 


We now consider methods for increasing the SNR, es¬ 
pecially for three-point correlations. The amount by which 
the SNR (3) must be enhaced is empirically estimated from 
the (problematic) example of * = 45°, whose reconstruction 
is more similar to a pristine image with i = 0° with a cor¬ 
responding MSE = 1.24 x 10 -3 . As shown in Figure O we 
have found empirically that SNR ( 3 ) needs to be increased 
by a factor of ~ 30 (corresponding to SNR ( 3 ) ~ 90) in order 
to see a noticeable improvement in the image reconstruction 
quality. 

One alternative to increase the SNR is to simply use 
much faster detectors, but even if much faster detectors ex¬ 
isted, their use would require the optical precision compa¬ 
rable to amplitude interferometers. We have also considered 
increasing the telescope area, but the telescope area cannot 
increase much more before individual telescopes would start 
to resolve the sou rce, resulting in a smaller correlation signal 
jRou et al.ll2013l ). 

Last we investigate the improvement of image recon¬ 
structions when using multiple spectral channels. For the 
fast-rotating star studied in this paper ( m v = 3, 50hrs), 
we have found that at least ~ 1000 spectral channels would 
need to be used in order to see a noticeable improvement 
in image reconstructions. When 1000 spectral channels are 
used with three-point correlations, we also find that recon¬ 
structed images are most similar to their corresponding pris¬ 
tine images, as shown in Table [f] In other words, Table [4] 
shows that the minimum MSE on the diagonal, unlike Ta- 
ble [3] Our lower limit of (~ 1000) required spectral channels 
is based on the assumption that the stellar image does not 
change much within the optical bandwidth, which is reason¬ 
able if we avoid observing over spectral lines. The required 
number of spectral channels may also be different for another 
telescope array configuration, e.g. one containing many re¬ 
dundant telescope pairs and triplets. 


6 DISCUSSION AND CONCLUSIONS 

Fast rotation is a common characteristic of hot, massive 
stars. It influences the structure and evolution during the 
whole lifetime of these stars, governing the loss of radiation, 
mass, and angular momentum. The circumstellar environ¬ 
ment of fast-rotators is thus intimately connected to the an¬ 
gular momentum in the central object. In the stellar interior, 
rotation induces circulations of mass, providing additional 
mixing of material that changes the stellar structure, life¬ 


time, production of chemical elements. The final stages of 
massive stars is also dependent on the previous evolution 
history, where rotation plays an important role. 

In the present work we demonstrate the capacity of SII 
to resolve the apparent intensity distribution of fast-rotating 
stars allowing to reconstruct high angular resolution images 
that directly reveal, independently of any adopted model, 
their position angle, apparent rotation flattening, and grav¬ 
ity darkening (polar to equatorial intensity distribution). 

The use of model fitting procedure can in addition 
be adopted to complement the reconstruction imaging ap¬ 
proach in order to estimate several stellar physical parame¬ 
ters and their uncertainties. 

The capacity of high angular resolution techniques to 
study fast-rotating stars was proven by recent results on 
image reconstruction and model fitting obtained from am¬ 
plitude interferometry in the last decade for the few bright¬ 
est (and thus closest) fast-rotating stars, which amount to 
half a dozen objects. 

The very high angular resolution attained by SII can 
take us one step forward allowing a more complete and sys¬ 
tematic study of stellar rotation effects on a much larger 
number of targets. For example, considering SII observa¬ 
tions at CTA with baselines up to 1000 m at wavelengths 
of ~ 500 nm allow to measure at least the three first vis¬ 
ibility lobes of stars having angular diameters as small as 
~ 0.5 mas. If the image reconstruction strategy presented 
here is adopted, we showed that for fast rotators brighter 
than m v = 3 a very large number (~ 1000) of spectral 
channels are needed in order for three-point correlations to 
be beneficial. However, improvements in the data analysis 
method may relax this requirement on the number of spec¬ 
tral channels. Even with these requirements, the number 
of possible targets for S II imaging studi es of rapid rotators 
given in the catalogue of lvan Belief (120121 ) is 35 (5 times more 
than the present number of observed fast-rotators). If one re¬ 
laxes the m v — 3 constraint then many more targets (nearly 
all the 354 stars in the van Belle’s catalogue, with roughly 
half this number per hemisphere) can be investigated by SII 
with simply two-point correlations (with the associated pos¬ 
sible complications shown in this work). Thus, compared to 
the performance of present amplitude interferometers, SII 
can potentially increase by a factor ~ 20 — 30 the number 
of observable targets for studies of fast-rotation across the 
H-R diagram. 

Finally, we note that although in this paper we con¬ 
centrated on hot, massive, fast-rotating stars, which is one 
of the most obvious science cases for SII with CTA, there 
are several other targets (e.g. binary stars, stars with photo- 
spheric spots from different origins) that are bright (m v < 3) 
and may have even stronger three-point correlation signals. 
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